Fermionization of two distinguishable fermions 
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In this work we study a system of two distinguishable fermions in a ID harmonic potential. This 
system has the exceptional property that there is an analytic solution for arbitrary values of the 
interparticle interaction. We tune the interaction strength via a magnetic offset field and compare the 
measured properties of the system to the theoretical prediction. At the point where the interaction 
strength diverges, the energy and square of the wave function for two distinguishable particles are 
the same as for a system of two identical fermions. This is referred to as fermionization. We 
have observed this phenomenon by directly comparing two distinguishable fermions with diverging 
interaction strength with two identical fermions in the same potential. We observe good agreement 
between experiment and theory. By adding one or more particles our system can be used as a 
quantum simulator for more complex few-body systems where no theoretical solution is available. 
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A powerful tool for solving complex quantum systems 
is to map their properties onto systems with simpler so- 
lutions. For interacting bosons in one dimension there is 
a one-to-one correspondence of the energy and the square 



of the wave function \tpixi 



to a system of iden- 



tical fermions [T|. The many-body properties of such ID 
systems have been studied in experiments with repul- 
sively interacting bosonic gases 0, Q . 
However, the essential property of a Tonks-Girardeau gas 

- namely the fermionization [l|, Ij] - is already present in a 
system of two interacting particles, regardless of the par- 
ticles being identical bosons or distinguishable fermions 
0. This two-particle problem is of significant interest 
because it is the main building block of all ID quan- 
tum systems with short-range interactions. It is also 
one of the few quantum mechanical systems for which 
an analytic solution exists. In contrast to measurements 
of bulk properties such as compressibility or collective 
oscillations [1, H, Q, we directly access the energy and 
the square of the wave function of the fundamental two- 
particle system. 

We realize such a two-particle system with tunable in- 
teraction using two fermionic ^Li atoms in the ground 
state of a potential created by an optical dipole trap and 
a magnetic field gradient (Fig. la, b). We can prepare 
this state with a fidelity of (93 ± 2)% |;7]. The energy 
of such two particles interacting via contact interaction 

- which is fully described by one parameter, the ID cou- 
pling strength g - can be analytically calculated for a 
harmonically trapped ID system The problem can 
be separated into center-of-mass and relative motion be- 
cause of the harmonic trapping potential and because 
the interaction term only depends on the relative dis- 
tance between the two particles. Then the solution can 
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FIG. 1: Trap setup and sketch of the performed experiment, 
(a) Our trap consists of an optical potential created by a tight 
focus of a laser beam and a magnetic field gradient, (b) De- 
terministic preparation of two fermions in the ground state 
of a potential well, (c) We measure the tunneling dynam- 
ics through a potential barrier for a repulsively interacting 
system of two distinguishable fermions for various interaction 
energies. The mean interaction energy per particle is indi- 
cated by the parameter U. These results are then compared 
with the tunneling dynamics of two non-interacting identical 
fermions in the same potential. 



be written as a product of the center-of-mass and the 
relative wave function. The latter is shown in Fig. 2a for 
different values of the coupling strength. For diverging 
coupling strength the square of the wave function of a 
system of two distinguishable fermions is the same as for 
two non-interacting identical fermions. This is the point 
where fermionization occurs. 

In our setup the particles are confined in a three- 
dimensional cigar-shaped potential with an aspect ratio 
of about 1:10, which can be harmonically approximated 
with trap frequencies of = 27r x (1.234 ± 0.012) kHz in 
longitudinal direction and ijJ± — 2tt x (11.88 ± 0.22) kHz 
in perpendicular direction. It has been shown in [9| that 
the energy of two interacting particles in the ground state 
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FIG. 2: Two particles in a ID harmonic potential, (a) Rela- 
tive wave function of two interacting fermions (blue) and two 
identical fermions (green) in a ID harmonic potential. For 
infinitely strong interaction (-1/(7|^|> 0) the probability to 
find the two distinguishable fermions at the same position 
vanishes. In this case the square of the total wave function of 
two distinguishable fermions is the same as for two identical 
fermions. (b) Relative kinetic energy. The blue and black 
curves show the energy of two interacting fermions in state 
I'll) depending on the coupling strength g\-[_i) given in units 
of ay = y^h/nLO\\ . The green line shows the energy of two 
identical fermions in state |tt)- The energy is plotted versus - 
l/g\^l) for a better comparison with the experimental results. 
The experimentally studied region is indicated by the dashed 
rectangle. 

of such apotential is well described by the ID solution 
given in [8|. Hence we treat our system in this ID frame- 
work. The combined optical and magnetic potential in 
one-dimensional form reads: 

VrMz)^pVo{l-Y:^-^J^)~^^rnB'z, (1) 

where Vq = ks 3.326 /iK is the initial depth of the optical 
potential, p is the optical trap depth in units of the initial 
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trap depth, zn — is the Rayleigh range of the optical 
trapping beam with minimal waist wq = 1.838 /im and 
wavelength A = 1064 nm, /i™ is the magnetic moment of 
the atoms and B' — 18.92 G/cm is the strength of the 
magnetic field gradient. The determination of the trap 
parameters is described in the supplementary material. 
The ID coupling constant g can be calculated from the 
3D scattering length a^o and depends strongly on the 
confining potential, which is characterized by the har- 
monic oscillator length a±_ = yJhfJujjZ [loj, where ft is 
the reduced Planck constant and — ^ the reduced 
mass of two ^Li atoms with mass m. The coupling con- 
stant is given by 



with C=C(^)=1.46... and C the Riemann zeta function. 
The value of g can be changed by tuning the 3D scat- 
tering length via a magnetic Feshbach resonance [U ■ 
When a^D approaches the extension of the confining har- 
monic oscillator potential a±^ , a confinement-induced res- 
onance (CIR) occurs for asD = a±/C 13]. Fig. 4b shows 
g\H) for two distinguishable atoms in the two lowest ®Li 
hyperfine states, |F= ■^,mp—--^) and |F= ■^,mp—-^) - la- 
beled It) and ID - as a function of the magnetic offset 
field fl4?|. For two identical fermions s-wave scattering is 
forbidden and thus g\^^^ = for all values of the mag- 
netic offset field. 

To determine the energy of the two-particle system in 
state Iti) we modify the trapping potential such that 
there is a potential barrier of fixed height through which 
the particles can tunnel out of the trap (see Fig. Ic and 
supplementary material) . In the presence of repulsive in- 
teractions the energy of the system is increased according 
to the blue curve in Fig. 2b. This decreases the effective 
height of the barrier and the particles tunnel faster. We 
allow the particles to tunnel out of the trap for different 
durations and record the number of particles remaining 
in the trap. By choosing an adequate barrier height we 
ensure that the time scale for tunneling is smaller than 
the lifetime of our samples in the ground state (about 
60s). Additionally obtaining meaningful tunneling time 
constants requires the timescale of the tunneling to be 
much larger than the inverse longitudinal trap frequen- 
cies of 0.7 ms. By averaging over many experimental re- 
alizations we obtain the expectation value of the particle 
number in the potential for different hold times (Fig. 3). 
By performing this measurement for various values of 
the coupling strength we can determine the dependence 
of the system's energy on g\^i). 

We find that for the observed range of interaction ener- 
gies - which are on the order of ftwy - only one particle 
leaves the potential even for long hold times. In a simple 
picture this can be explained as follows: If one particle 
tunnels through the barrier the interaction energy is re- 
leased as kinetic energy, which leaves the other particle 
in the unperturbed ground state of the potential. This 
state has a tunneling time scale much larger than the du- 
ration of the experiment. Thus we can fit exponentials 
of the form N{t) — iVtunncie~~ -I- A^remain to the mean 
particle number to deduce the tunneling time constant 
T for different magnetic fields. The mean numbers of 
tunneled (iVtunnoi) and remaining particles (A^^remain) are 
expected to be unity. However, due to the finite prepara- 
tion fidelity they are slightly lower. In Fig. 4 we show the 
determined tunneling time constants of a system of two 
interacting fermions for different interaction energies as 
a function of the magnetic field. We observe a decrease 
in the tunneling time constant over two orders of mag- 
nitude for increasing magnetic field due to the gain in 
interaction energy caused by the CIR. 
For a direct comparison of the properties of the two in- 
teracting distinguishable fermions with those of two iden- 
tical fermions we perform the same measurement with 



3 



distinguishable fermions, 750G (g^^ > 0) 
distinguishable fermions, 782G (|g^J » 0) 
distinguishable fermions, 900G (g^^ < 0) . 
identical fermions, 782G 
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FIG. 3: Mean number of particles remaining in the poten- 
tial well. After modifying the initial potential the particles 
can tunnel through a barrier of fixed height for a certain hold 
time. Subsequently, tunneling is switched off and the mean 
particle number left in the potential is recorded by averaging 
over many experimental realizations. Exponential fits to the 
data (solid lines) allow to extract the tunneling time constants 
of two interacting distinguishable fermions for different inter- 
action strengths (blue) and of two identical fermions (green). 
Each data point is the average of about 70 measurements ex- 
cept for the first and the last data point in each series (about 
230 realizations). The errors are the SEM. 



two fermions in state |f|) in the same potential (Fig. Ic). 
The results of these reference measurements are shown in 
Fig. 3 and Fig. 4 (green points). As the identical fermions 
are non-interacting we find no dependence of the tunnel- 
ing time constant on the magnetic field in this measure- 
ment. 

Comparing the results of the two systems we find that 
the tunneling time constant for the interacting system de- 
creases monotonically with increasing magnetic field and 
crosses the magnetic field independent tunneling time 
constant of the two identical fermions. Thus there is one 
magnetic field value where the tunneling time constants 
of both systems are equal. At this point both systems 
must have the same energy. For a ID system with given 
energy there is only one unique solution for the square 
of the wave function. Therefore, right at the observed 
crossing point of the tunneling time constants the en- 
ergy and the square of the wave function |V'(2:i,Z2)P of 
the two interacting distinguishable fermions and the two 
non-interacting identical fermions must be equal. Hence, 
exactly at this crossing point the system of two distin- 
guishable fermions is fermionized. As predicted by the- 
ory Q we find the position of the fermionization at 
the magnetic field value where g\fi) diverges due to the 
confinement induced resonance. 

For magnetic field values below the CIR we have realized 
the two-particle limit of a Tonks-Girardeau gas. Above 
the CIR we have created a Super- Tonks state consist- 
ing of two particles. The Super- Tonks state is a strongly 
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FIG. 4: (a) Tunneling time constants for different values of the 
ID coupling strength. The tunneling time constant of two re- 
pulsively interacting distinguishable fermions (blue curve) de- 
creases by two orders of magnitude with increasing magnetic 
field. This is attributed to the gain in interaction energy when 
ramping across the CIR. The tunneling time constant of two 
non-interacting identical fermions (green line) remains unaf- 
fected by the magnetic field within our experimental accuracy. 
At the magnetic field value where both curves cross we iden- 
tify the fermionization of two distinguishable fermions. The 
errors are the statistical errors of the fits shown in Fig. 3. The 
blue line is a guide to the eye. (b) One-dimensional coupling 
constant g\^i) with a CIR at (783.4 ± 0.4) G. For the calcu- 
lation we used the perpendicular harmonic oscillator length 
a± = i/fi/ imjj± of the modified potential. 



correlated metastable state above the attractive ground 
state branch (see Fig. 2b). In a system with particle num- 
bers > 3 inelastic three-body collisions lead to a fast 
decay of the metastable Super- Tonks-Girardeau gas [s!]. 
In contrast, our two-particle Super- Tonks state is stable 
against collisional losses since there is no third particle 
available to undergo an inelastic three-body event. 
To determine the energy of the two interacting fermions 
from the measured tunneling time constants we use a 
WKB calculation (see supplementary material) . This re- 
quires knowledge of the potential shape. The parameters 
of the optical potential are determined by precise mea- 
surements of the level spacings in the potential. The 
final parameter to determine the barrier height is fixed 
by the measured tunneling time constant of two identi- 
cal fermions (see supplementary material). The energies 
obtained from the tunneling time constants of two dis- 
tinguishable fermions are shown in Fig. 5. 

We compare these energies to the analytic theory for 
a harmonic potential [8] (see Fig. 2). This theory needs 
two input parameters, the coupling strength and the level 
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FIG. 5: Interaction energy of two fermions for different inter- 
action regimes. By using a WKB based calculation we can 
determine the energy of two distinguishable fermions at dif- 
ferent interaction strengths (blue points) from the tunneling 
time constants presented in Fig. 4a. The blue curve shows 
the expected energy shift for a harmonically trapped system 
(dashed rectangle in Fig. 2). 

spacing. For the coupling strength we use g|-f-j^) of our sys- 
tem shown in Fig. 4b. For the level spacing we use the 
energy difference fi.a;|| ^aic = Eq — Ei = 2^1% x 743 Hz be- 
tween the ground and first excited state of the potential 
which we calculate using the WKB method. With this 
approximation the energy obtained from the tunneling 
measurements and the energy obtained from the ana- 
lytic theory Q are the same at the CIR. For the Tonks 
regime we find excellent agreement of the experimentally 
determined energy with the theoretical prediction for a 
harmonic trap. Above the CIR the harmonic theory is 
not applicable because the second excited state is not 
bound in our potential. Additionally, we expect devi- 



ations for larger energies due to the limited validity of 
the WKB approximation for energies close to the con- 
tinuum threshold. A more precise description could be 
achieved by adapting the theory described in ^ to our 
non-harmonic potential and by usirig a more accurate 
theory for the tunneling process [H, . 

In summary, we have measured the interaction energy 
of two distinguishable fermions as a function of the inter- 
action strength and identified the point of fermionization. 
The good agreement between our results and theoretical 
predictions shows that our experiment has the capability 
to simulate strongly correlated few-body quantum sys- 
tems. Using the experimental methods established in 
this work it is straightforward to extend our studies to 
more complex systems. Simply addin g a third particle 
either in one of the present spin states [13| or a different 
spin state [H, [l^ allows us to study a highly non-trivial 
system where no analytical solution exists. In a few-body 
system with defined particle number and attractive inter- 
action we could investigate pairing phenomena and thus 
work towards studying superfluidity in finite systems. 
This has already been investigated in the context of nu- 
clear physics [2Q|] . By dynamically changing the shape of 
the trapping potential we could simulate a vast amount 
of different time-dependent quantum systems. A feasi- 
ble experiment would be to periodically modulate the 
strength of the magnetic field gradient. This would allow 
us to study ionization-like excitations in the strong-field 
regime |2ll | which have been studied in ultrafast physics 
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Supplementary material 

Tunneling potential. Our initial samples, which 
we prepare as described in Q, consist of two non- 
interacting particles in the ground state of our trap. 
At the end of the preparation the optical trap depth 
is p — 0.795 (relative error of p: 1.3 x 10"'^) and 
thus there are 4 bound states in the potential. At 
this trap depth tunneling from the two lowest levels 
is completely suppressed. For the measurement of the 
interaction-induced tunneling we ramp to a trap depth 
oi p = 0.6875 with a ramp speed of dp/dt — 0.043 ms~^ 
where tunneling occurs on experimentally accessible 
time scales. In Q we have estimated the probability 
of exciting particles when performing a single ramp 
of the optical potential at this speed to be consistent 
with (3 ± 1)%. To deduce the tunneling time constant 
we fit the measured particle number with a X"Square 
minimization algorithm (Levenberg-Marquardt) con- 
sidering the statistical error of the mean particle number. 

Magnetic offset field. To tune the strength of the 
interaction we apply magnetic offset fields ranging from 
523 G to 900 G. At two distinct magnetic field values 
close to the CIR (783.2 G and 788.9 G at p = 0.795) 
we observe molecule formation resonances with a width 
of 0.4 G (FWHM). We attribute these resonances to a 
coupling between the relative motion of the two atoms 
and the center-of mass motion of a molecular bound 
state due to the anharmonicity of the trapping potential 
[23I [23 ] . We prevent formation of molecules by ramping 
across those resonances with sufficiently high speed (20 
G/ms). Thus we can neglect these resonances for our 
experiments. 

WKB approximation. To determine the energy 
of the bound states of the potential with the WKB 
(Wentzel-Kramers-Brillouin) approximation one has to 
solve the following implicit equation for the energy E: 

■\ f ' y^2^^U^^^V{^ dz = n+l 

with n a positive integer number and Za < -Zfa < -^c 
the solutions of V{z) — E = 0. The tunneling time 
constant r of a particle with energy E is calcu- 
lated by T = / /T with the transmission coefficient 



T = exp{-2 Q ^^{V{z)- E)dz) and the knock- 
frequency I = For the energy range studied in 
the experiment we observe that only one of the two 
particles tunnels through the barrier. This particle 
has an energy identical to the relative kinetic energy 
of the two-particle system. Hence we can determine 
the relative kinetic energy of the two-particle system 
by extracting the energy of the tunneled particle from 
the tunneling time constant Te^^p. To map r^xp onto 
energies we calculate tunneling time constants for a set 
of energies. By matching these calculated tunneling 
time constants to the measured ones we determine the 
relative kinetic energy of the two-particle systems. A 
more detailed description of this method will be given in 

Determination of the potential shape. The 

potential consists of a cigar-shaped cylindrically sym- 
metric optical potential Vopt created by a tightly focused 
laser beam and a linear magnetic potential Vmag- The 
longitudinal part of the potential is given by Vr=o{p^ z) = 

Vopt{p,z) + Vmag{z) = pVq (1 - (i + (^/^^,)2)) ) " Mm B'z 

where Vq is the initial depth at the center of the optical 
dipole trap, p the optical trap depth as a fraction of 
the initial depth, z^ = vrwg/A the Rayleigh range, 
A = 1064 nm the wavelength of the trapping light, B' the 
magnetic field gradient and the magnetic moment of 
the atoms. 

The optical potential is created by a tightly focused 
laser beam with a power of Pq = (265 ± 27) //W. We 
assume that the beam has a Gaussian shape, neglect- 
ing any aberrations that might be introduced by the 
optical setup. To deduce the shape of the potential 
we directly measure the structure of the energy levels 
in the optical potential by periodically modulating the 
power or position of the trapping beam. For this we use 
non-interacting atoms and thus the observed excitation 
frequencies correspond to the energy differences between 
the different trap levels. At an optical trap depth 
of p = 1 we find level spacings of huj\\g__-^ = 27rfi, x 
(1.486 ±0.011) kHz, ^ 2nh x (2.985 ± 0.010) kHz 

and fta;||,_^ = 27rft x (2.897 ± 0.020) kHz for the lowest 
trap states in longitudinal direction where the indices 
label the energy levels in the trap. To obtain the waist 
we vary wq and Pq within its error crp„ to minimize 

S(i,j) = (0,l)(0,2)(2,4) ~ [{^optj - Eopti) - 

Wi-j 

with Eopti j being the energy of the calculated bound 
states of the varied optical potential. The bound state 
energies are calculated using a WKB calculation. We 
find a minimum deviation for wq = 1.838/im and 
Po = 291.5 /iW resulting in an initial depth at the 
center of the optical trap of Vq = fc^ x 3.326 /iK. For 
these parameters two of the three frequencies match the 
calculated bound states within their errors and all three 
match within 2a. After fixing the parameters for the 
optical potential we have to determine the strength of 
the linear magnetic potential which is created by a mag- 
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netic field gradient B' . From a levitation measurement 
we obtain B' = (18.9 ± 0.2) G/cm. To obtain a more 
accurate parametrization of the potential barrier we use 
the tunneling measurement of two identical fermions at 
a trap depth of p = 0.6875 as a calibration. We perform 
a WKB calculation of the tunneling rate and modify 
the value of B' in the calculation such that the resulting 
tunneling time constant agrees with the experimentally 
observed tunneling time constant of (74.1 ± 2.7) ms. 
From this we obtain B' = 18.92 G/cm. Using this 
magnetic field gradient the calculation is consistent 
with the deterministic preparation [7| of (2, 4, 6, 8, 10) 
particles at p=(0.6575, 0.7025, 0.7475, 0.7863, 0.8200) 
in the experiment. 

We measure the perpendicular trap frequency by exciting 
non-interacting particles in the pure optical potential 



Vz=o {r, p) = pVq e~ . In perpendicular direction we 
observe two resonances at 2ijJx — 2tt x 29.07kHz and 
2ijjy = 2tt X 28.26 kHz for p = I which we attribute to 
a slight ellipticity of the trap. Following [H [13, [11] we 
only consider one CIR which we calculate from the mean 
perpendicular frequency of 27r x (14.33 ± 0.26) kHz. 
At p = 0.6875, where we perform the tunneling mea- 
surements, the trap frequency of the perpendicular 
confinement has to be scaled by and is given by 
a;_L = 27r X (11.88 ± 0.22) kHz. At this depth the mean 
longitudinal trap frequency of the pure optical potential 
is given by uj\\ = 2-k x (1.234 ± 0.012) kHz which is 
calculated from the excitation frequencies '^||o_i and 

^ll0-2- 



